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In this paper, instabilities of the isotropic metallic phase in iron pnictides are investigated. The 
relevant quartic fermionic interaction terms in the model are identified using phase space arguments. 
Using the functional integral formalism, a Hubbard-Stratonovich transformation is used to decouple 
these quartic terms. This procedure introduces several bosonic fields which describe the low-energy 
collective modes of the system. By studying the behavior of these collective modes, a possible 
instability is found in the forward scattering channel of the isotropic phase driven by magnetic 
fluctuations. Using mean field analysis, we obtain a static and homogeneous ground state. This 
ground state is metallic, but the electron Fermi pockets are distorted unequally at different pockets in 
momentum space. This results in a desirable nematic ordering which breaks the lattice d symmetry 
but preserves translational symmetry and may explain several experimental observations. 
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INTRODUCTION 



A new unconventional superconducting state 
has recently been discovered in doped iron pnic- 
tide (FeAs) materials }^- Although FeAs is sim- 
ilar to the cuprates in some respects, i.e. it is 
a layered and exhibits magnetic ordering in its 
parent compound, some of its electronic proper- 
ties are fundamentally different. In contrast to 
cuprates, the electronic structure of FeAs is ef- 
fectively described by a multi-band model4i^— In 
addition, FeAs is believed to be less correlated 
than the cuprates. These observations suggest 
that their mechanism of superconductivity may 
be different. In FeAs, the Fermi surface (FS) 
breaks down into several pockets. Due to this 
feature, scattering processes at low energy can 
give rise to non-trivial many-body physics. This 
may explain the origin of superconductivity in 
FeAs with only repulsive interactions 4 

In addition to the superconducting state, an- 
other intriguing property in FeAs has recently 
been observed. The electronic properties ex- 
hibit a directional preference, ~"iS. regardless of 
the system's intrinsic 4-fold crystal rotational 
symmetry (C4 symmetry) . This anisotropy may 
be related to nematicity due to the interac- 
tions between electronsjii^ The new phase is 
metallic and associates with magnetic order. 
This observation has led to speculations about 



the relationship between magnetic fluctuations 
and nematic ordering. Some theoretical stud- 
ies have shown that magnetic fluctuations can 
play an important role in giving rise to ne- 
matic orderiii"— The functional renormaliza- 
tion group (FRG) point of view2ii2^ provides an 
interesting perspective to understand the stabil- 
ity of the metallic phase in FeAs. Because of the 
intricate geometrical structure of the FS, the 
coupling parameters in the forward scattering 
channels are relevant under RG transformation 
instead of just being marginal. This suggests 
that some non-trivial scattering processes may 
occur in the direct channels which can poten- 
tially break the stability of the FS. 

Based on these observations, the central goal 
of this paper is to investigate how magnetic fluc- 
tuations can affect the direct channels in FeAs, 
and yield instabilities. We argue that nematic 
order is more favorable energetically than the 
isotropic metallic phase. This paper is orga- 
nized as the following. In Sec. |lTl we will dis- 
cuss the effective model that will be used in the 
study. Sec. IIIII will show how the magnetic 
fluctuations affect the density fluctuations and 
give rise to instabilities in the isotropic metallic 
phase. Sec. IIVI the nematic order and its mag- 
netic fluctuations will be discussed at the mean 
field level. In Sec. |Vl we will briefly discuss 
these results in the context of experiments and 
other theoretical studies. 



II. MODEL HAMILTONIAN 

We extend the model of Ref . [20 by including 
the forward scattering (density-density) interac- 
tions within intra- and inter- FS pockets. The 
free part of the Hamiltonian is given by 



^o = E^£^L 



ks^a.ks, 



(1) 



where s is the spin, and a represents the lo- 
cation of the FS pockets, a = F represents 
the hole FS pocket with its center located at 
(0,0), and a ~ X,Y represent the electron FS 
pockets with centers located at Qi — {it, 0) 
and Q2 = (0, tt) respectively (see Fig. [T]). 
The energy dispersions are e^ = co ^ t /^, 
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= —en + 7^ — h Tp^ w, and e?f , r^ 



_ kt . kt 

^*^n + 2^ + 2^ ~ M' where m, m^, and ruy 
are the band masses, eo is the offset energy, and 
^ is the chemical potential. For simplicity, we 
shift the momentum of the X (Y) pocket to Qi 



(Q2) i.e. e^+Q^ -^ e? (el^r^^ -^ el). Fur- 
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thermore, summing over all spins is implicitly 
assumed throughout the paper. 
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FIG. 1. A schematic diagram showing the FS struc- 
ture of a typical FeAS in the unfolded Brillioun 
Zone: The bold curves represent the FS of the 
model in the isotropic metallic phase, the dot- 
ted curves represent the unequally renormalized FS 
pockets in the nematic phase with (A^) > 0. 

In the low temperature limit, the dominant 
scattering processes only take place in the vicin- 
ity of the FS. Moreover, because of momentum 



conservation, all scattering processes are strictly 
constrained to the available phase space for two- 
body scattering processes. We are interested 
in the processes with small momentum transfer 
(or up to a nesting vector in exchange chan- 
nels), as these processes are likely to develop 
collective modesJS, Therefore, keeping only di- 
rect, exchange, and BCS channels in the model 
is sufficient, since these dominate over the oth- 
ers in phase space. 

Although most of the momentum channels 
are being discarded, there is still ambiguity 
which of these three special channels are rel- 
evant to the problem. We can further elim- 
inate some of them by using the following 
intuitive arguments. First, for electron-hole 
(e-h) interactions, because of the FS pocket 
nesting spin-density-wave (SDW) fluctuations 
are the most important collective mode, espe- 
cially as the system temperature approaches 
the magnetic ordering phase transition temper- 
ature (T/v). Therefore, the exchange channel in 
the e-h interactions is kept. Normally, the ex- 
change channel of this interaction breaks into 
spin triplet and singlet channels, which corre- 
spond to SDW and charge-density-wave (CDW) 
respectively. However, the CDW fluctuations 
are insignificant and can be discarded, because 
their characteristic energy scale is much lower 
than for the SDW. By using similar arguments, 
all the direct and BCS channels in the e-h in- 
teractions can be discarded. For the electron- 
electron (e-e) interactions, since we wish to in- 
vestigate the instabilities of the isotropic para- 
magnetic phase (metallic), the direct channels 
should be kept. Following the same reasoning 
for the e-h interactions case, the exchange and 
BCS channels are irrelevant in e-e interactions. 
Therefore, the effective interacting Hamiltonian 
becomes 



^int — 



^spin 



/ . / . ^".q ' ^".-q 



a=X,Y q 



(1) 



a=X.Y q 

q 
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where s, 



a,q 



= E, 



J 



the Pauli 

J 



k^r,k+q,s^ 

matrix vector, 



__' Casks', cr IS 
and pa,q = 
I]kcL,k+q,sCa,ks- The first term of Eq. ^ 
combined with "Ho is the original model given in 



to check whether u\ ~ uq is positive or 
negative. Moreover, we will see later that this 
is also the key parameter determining whether 
instabilities can emerge in the direct channels. 



Ref. |20|. Uspin is the coupling constant for elec- For now, we assume 



,W 



Uq is negative 



(1) 
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tron and hole pocket interactions in the triplet 
channels. The last two terms are newly intro- 
duced, describing the density-density interac- 
tions between the X and Y FS pockets, u 
is the coupling constant for the intra-electron 
pocket interactions, and wg is the coupling con- 
stant for inter-electron pocket interactions. The 
notation of the coupling constants in ^ are the 
same as in Ref. [23, and angular dependence has 
been ignored. Because of the C4 symmetry in 
the system, the electron intra pocket interac- 
tions (second term in ©) have the same cou- 
pling constants. The above interaction terms 
only capture dominant scattering processes of 
the electrons in the vicinity of the FS. The full 
interaction terms include arbitrary momentum 
transfer contributions as well.— 



and set ue 



,(1) 



and ue 



,(1) 



(positive) . Then we introduce a set of auxiliary 
bosonic fields (j)q = (A+, A", M^, M^, where 
q — (ii/„,q) and ^'„ is the bosonic Matsubara 
frequency. 

By inserting the 'fat' unity (jAip into the par- 
tition function integrand {Z = J X'[^, ^t^Jle""^), 
the fermonic quartic terms in the action are 
canceled out by shifting the auxiliary fields 
(jA2l - [X4|) . By introducing the Nambu spinor 

*fe = i^kv'^kv^lt''^kv^lv'^ki)^ the new ac- 
tion can be compactly expressed as 



5 = E ^'^^kk'^^'k' + -^ E M^ M^, 



kk' 



E 



tA;a^ 






(3) 



III. INSTABILITIES OF THE 
ISOTROPIC PHASE 



where 



To investigate whether the model in Sec. |TI] 
can develop any instability, one can use rig- 
orous diagrammatic methods. However, we 
use a simple alternative approach to tackle the 
problem. This approach is based on the func- 
tional integration formalism, introducing aux- 
iliary bosonic fields to decouple the fermonic 
quartic terms (Hubbard- Stratonovich transfor- 
mation). By integrating out the fermonic de- 
grees of freedom, one can obtain the action of 
these bosonic fields which describes the collec- 
tive modes of the system at low energy. In the 
following, we demonstrate this procedure and 
investigate the stabilities of the isotropic metal- 
lic phase. 



A. Hubbard-Stratonovich Transformation 

In order to perform the Hubbard- 
Stratonovich transformation, it is important 



Skk' =5o,fcfc' +Vfcfc'. 



So,kk' is given by 
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6kk' 



and Vkk' is given by 



^iA 



k-k' 



^k-k' 



Mf_,, 



^^t-k' 



A 



MLk' 



k-k' 



(4) 



(5) 



Mtk' 

ML,, 





(6) 
where G^\ — — iw„ -|~ eg, and a ~ X,Y, T. 

A~ must change sign under exchanging x and 
y axis to preserve the original C4 symmetry in 
([3]) because the boson-fermion coupling terms. 
Integrating out the fermionic fields, we obtain 



the action for the bosonic fields, 



a.,q 



^Ef;^A^At, + i-A-A=, 



2^^ \u 
1 



(7) 



Note that M" ~ M"n, where n is some arbi- 
trary direction of the spin. Since A^ change 
sign under C4 rotation, any non-zero ground 
state expectation value of A^ implies sponta- 
neous C4 rotational symmetry breaking. There- 
fore, we can identify A~ as the nematic order 
parameter, and Eq. ([7]) describes how the mag- 
netic and nematic collective modes interact with 
each other. 



B. Isotropic metallic phase 



The extremum of the action is given by 



S[^] 



(8) 



q = (<t>q) 



Solving this equation involves the inversion of 
Q~^ . However, without knowing the functional 
form of (t)q, this is a formidable task, since Q~^ is 
generally not diagonal in momentum-frequency 
space. Therefore, in order to proceed, we need 
to guess a solution which simultaneously sat- 
isfies Y.qQkq[Q^%k' = I3x34fc' and Eq. ©. 
Since we want to investigate the stability of the 
metallic phase, one obvious and desirable choice 
is (M^) — (M^) — (paramagnetic), and 
(A^) = (isotropic). In general (A+) can be a 
nonzero constant. For instance^ the FS pock- 
ets can shrink or expand according to its non- 
zero expectation values. However, the hole and 
electron FS pockets must be renormalized, such 
that no net charge is introduced into the sys- 
tem (Luttinger theorem))^ Modifying the FS 
pockets in this manner does not influence our 
result qualitatively, since angular dependence in 
the coupling constants is ignored. Therefore, we 
assume (A+) = 0. 

Using the saddle point approximation, we 
expand Eq. ([7]) around the extremum with 



{(j)q) = and include fluctuations beyond the 
mean field level. Thus, we obtain 

trlng"^ -trln[^o-i(l-fgoV)] 

= trlng-i-5:i-tr(eoV)^", ^^^ 



and 



'(5A 



k~k' 



5Mtk- 



Vkk' = I -<5A,~_,, 5Ml_^, I , (10) 
^5M^_^, 5Ml_,, 

where 5M^^^, and 5A^_^, arc the fiuctuations 
of the collective modes near {(pq) = 0. This 
yields the effective action 



5, 



eff 



+ A5^(^A-)2(<5A/4")^ 



(11) 



where H, = Y.kGa±Ga,k+q-, Xa,g == :^ + 

^X^fc ^r,A;Ga,fe+g, and A = X^fc ^r,fe(Ga,fe) , 
and X = (r, x) represent the imaginary time and 
real space coordinate respectively. Note that 
the frequency and momentum dependence in A 
can be discarded because of power counting^^ as 
only local interactions are considered. ^^ More- 
over, the other terms beyond Gaussian only 
yield higher order corrections for the action and 
are not important to the problem (before the 
magnetic order sets in), except in the last term 
of (ITlT) . which accounts for the interactions be- 
tween magnetic and nematic order fluctuations. 
The physical meaning of this term is that mag- 
netic fiuctuations associated with the fiuctua- 
tions of the X and Y FS pockets unequally can 
lower the total energy, since this term is neg- 
ative (A is negative). If the magnetic fiuctua- 
tions are small, these contributions are negli- 
gible. When the magnetic fluctuations become 
large in the vicinity of the magnetic ordering 
transition point, the stability of the (A^) = 
isotropic phase can break down. 



Instabilities due to magnetic 
fluctuations 



To find the instabilities, we can integrate out 
5M". For small q, Xa}q = ro + jWn] + fa,q, 
where ro = Xa,q=o: 7 i^ the Landau damping 
coefficient, and fa,q = ql{l±ri)+ql{l^ri)+rizq'^^ 
is a generalized anisotropic function with — 1 < 
rj < 1, and the upper (lower) sign refers to 
a — X (a — Y)i^ In this paper, we are not in- 
terested in the out-of-plane anisotropy. There- 
fore, we set rjz to zero. We thus can obtain the 
effective action for i5A~, 






(12) 



4- tr ln[ro + -fdr + faA + H^A-f], 

where q are momentum operators. The trace 
in the above equation can be evaluated in 
momentum-frequency space, yielding 



a,g 



ln[ro + -flunl + fa.q + A(5Ag ^A_ J. (13) 



Now, expanding the logarithmic function with 
respect to (5A~)^, we can further simplify it by 
keeping only the leading term. 



•S'eff ' 



q 



Uq 



,,(1) 



A 



ro+l\l^n\ + fa 



2n„ 



5A: 



(14) 



As the temperature approaches Tat from above, 
the first two terms in the coefficient of this 
Gaussian term remain finite and positive, but 
the last term is negative and diverges. There- 
fore, the coefficient approaches zero and even- 
tually turns into the negative regime. This sig- 
nals a new instability, setting in prior to the 
divergence of the magnetic susceptibility. This 
instability occurs in the direct channels and it 
is similar to Pomeranchuk instabilities.^'* This 
implies a deformation of the FS in the ground 
state fi^ Also note that since 11^ > 0, this insta- 
bility cannot arise without the aid of magnetic 
fluctuations, and nesting between electron and 
hole pockets is essential in this scenario. 



IV. NEMATIC ORDERED PHASE 

So far, we have shown that the isotropic 
metallic phase can be unstable near Tjv due 
to magnetic fluctuations. This suggests that a 
new phase occurs between the isotropic metallic 
and the magnetic phase. In this section, we will 
identify this as nematic ordering. The effective 
action for the nematic order parameter will be 
derived within mean field theory. 



A. Mean field theory 

To obtain the Ginzburg-Landau action, one 
cannot simply replace SA^ by its order param- 
eter in (jlll) . This method works only when the 
order parameter is small and all the terms be- 
yond the Gaussian approximation are positive. 
Although the order parameter is small in our 
case, the quartic term (magnetic-nematic fluc- 
tuation coupling) is negative. Any truncation of 
the series expansion beyond the quartic terms 
in dni) will potentially lead to an ill defined (di- 
verging) functional integral. Therefore, instead 
of expanding around {(f>q) = in the (/)-field 
manifold, we only expand A+, M-^ and M^ 
around zero and keep A~ in an arbitrary func- 
tional form. Thus, Qq^ is replaced by 5'q , 



S'o.kk' = ^1 



0,kk' 








-a: 



, (15) 



and V is replaced by V", 



V'kk' = 



M'ltk' SM^-k' 



SM^~k' 

JiVC,, I . (16) 







Again, the fluctuations of A+ are ignored. 
Thus, 



SeS 



lE 



trlnC/'g 



-A^A" - 

- trln[l + e'o V] 



-Y^M-M-^ 



(17) 
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So far, the above equations have not assumed 
any functional form for A~. In order to pro- 
ceed, we assume that the stable ground state en- 
ergetically favors homogeneous configurations 
in space (preserving translational symmetry) 
and time (static). Then, we can invert G'q ex- 
actly and expand ([TT]) around (A+) = (M^) = 




•-•cff — "-"n 



where 



[A- 






SM^X. 



"~.-iJA-f" 



(18) 



Z Ur\ 



EMGx:fc + A-][G^,i-A- 



(19) 



and x^}^ = 2/u3pin + 2^^ l/[GFi(G-i,+^ ± 
A^)], where the upper (lower) sign refers to 
a ^ X {a = Y). The effective action in (fTSl) 
is our main result that describes the low-energy 
magnetic collective modes in the nematic phase. 
'5'ncm is the Ginzburg-Landau action of the ne- 
matic order parameter, and the second term 
describes the magnetic fluctuations. The origi- 
nal magnetic and nematic fluctuations coupling 
term in (jlip is implicitly contained in the sec- 
ond term of Eq. (fT8|). One can quickly check 
that in this phase the magnetic fluctuations are 
unequal in the X and Y direction. Namely, the 
susceptibility (the correlation of magnetic fluc- 
tuations) is {SM^SM^,) ^ Xa,k-k', which is A~ 
dependent. 

Furthermore, at this point the system still 
possesses a SDW instability, since Xa\ ^P" 
proaches zero at a finite temperature Tjy. But 
T/v niay be different from the prediction in 
isotropic phase, because non-zero values of A~ 
modify the nesting condition between the elec- 
tron and hole FS pocket. This modification lifts 
the degeneracies between stripe orders. 

Note that the magnetic fluctuations peak 
around Qi and Q2 or q = at low energies. 
Therefore, it is natural to impose a cutoff fre- 
quency Wc for SM" This frequency is also the 
characteristic frequency that describes the typi- 
cal wavelength of the magnetic fluctuations. To 



make a connection with the phenomenology of 
magnetic susceptibilities, we can constrain our- 
self to the fluctuating modes below this cutoff. 
Hence again, by expanding Xa^q around q = 0, 

we obtain Xa!q = Xa!q=0 + ^Wn\ + fa,q, but 7 

and fa,q can be different in the isotropic phase. 
If the FS renormalization is not too drastic in 
nematic phase, we can further expand Xa q=o 
around the isotropic phase. Thus, we have 
Xa q=o = ''o + kA^ , where k is some constant 
that depends on the microscopic details of the 
material. This result coincides with Ref. 123 . 

The mean field equation of the nematic or- 
der parameter can be obtained by varying S'nom 
with respect to A~, 



2A- 



E 



1 



1 



G- 



Gl 



A- 



'X,k ' " ^Y,k 

(20) 
This equation can be straightforwardly be 
solved numerically, giving the location of the 
minimum of 5nom- One can also quickly check 
that A^ = satisfies the mean field equation 
trivially. 



B. Nematic ordering phase transition 

To search for the non-zero of A~, instead of 
evaluating the mean field equation, we plot the 
value of 5nom[A^] versus A^ at a given temper- 
ature numerically. This demonstrate how S'nom 
can develop a minimum at non-zero A^ , as tem- 
perature is varied. In order to do this, we eval- 
uate the Matsubara sum in (fT9|) . 



Sn 



,[A-] = iA(A-)^ 
2 u„ 



k 



Infl-^e 



-0i< 



)][l + e- 



-PieZ^A-)^ 



(21) 



As shown in Sec. IIII[ the isotropic state is un- 
stable only when the magnetic fluctuations be- 
come important. Therefore, it is natural to take 
LUc as the cutoff energy for this effective model. 

Then,Ek^r^(27F'^^/o''^^^- 



The input parameters are 62 — 0.2eo, /i ~ 
0.05eoj^^ where eo — 0.2 eVj and the mag- 
netic fluctuations cutoff energy Wc = 3.4 meVi^ 
In addition, we set a — (l/M„)(27r)^(eo/m-), 
which is a dimensionless parameter controlling 
the transition temperature in the mean field ap- 
proximation. For demonstration purposes, we 
set Un = 0.1 and use the hole Fermi pocket mo- 
mentum from Ref. lS. Thus we have a ~ 1.5. 
Using these input parameters, we obtain the 
plot shownthe plot is in Fig. [5J 
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FIG. 2. The nematic order phase transition oc- 
curs at some finite temperature where S'nem = 
S'ncm(27r)V(meg). 



According to Fig. [21 the transition tem- 
perature from the mean field approximation is 
/3eo ~ 12 or T ~ 193 K and A" ~ 0.4eo. These 
values are clearly overestimates compared to ex- 
periments. This is because mean field theory 
does not taking account the correlations be- 
tween fiuctuations. This deviation is even more 
pronounced because in the low effective dimen- 
sionality. However, the mean field result at least 
estimates the correct order of magnitude. 

So far, magnetism does not appear to play 
any direct role in the nematic phase. However, 
the order parameter depends on the cutoff Wc 
(this is similar to the mean field analysis in 
BCS theory). This determines the position of 
the minimum, where by a smaller cutoff implies 
smaller A~ or lower transition temperature and 
vice versa. To further investigate the relation- 
ship between magnetic fiuctuations and nematic 
order calls for more sophisticated methods be- 
yond mean field theory. 



V. DISCUSSION 



A. The control parameter 



For the nematic order considered in this pa- 
per, u„ = ug — M4 is the key parameter. It not 
only determines whether the system can sustain 
nematic ordering, but also determines the tran- 
sition temperature to the nematic order from 
the isotropic phase. These coupling constants 
depend on microscopic details, and are difficult 
to obtain in general. Nevertheless, this could 
explain the fact that nematic ordering is not 
found in every FeAs material. 

Why Un plays such a crucial role in this sce- 
nario can be understood by the following in- 
tuitive picture, u^ and U4 are related to the 
coupling strength between two quasi-particles, 
depending on whether both reside in different 
electron FS pockets or the same electron FS 
pocket respectively. If u^ > u\ means that 
two quasi-particles tend to occupy the same 
electron pocket rather than different pockets 
because this yields a smaller potential energy. 
Therefore, instabilities of the original isotropic 
ground state may arise, whenever this scatter- 
ing processes surpass the other channels. 



B. Comparison with other work 

The nematic order considered in this study 
also yields an anisotropic magnetic fluctua- 
tion, which can explain results from recent 
neutron scattering experiment 1^ It also im- 
plies anisotropic resistivity in the nematic 
phase. i^^i^ 

In earlier theoretical work, it was shown ex- 
plicitly that nematic order can be developed 
solely from the effect of magnetic fluctuations, 
and yields an anisotropic magnetic fluctuat- 
ing state. In the phenomcnological approach 
based on a Ji — J2 modelji^r— or a preemp- 
tive nematic order from a microscopic itinerant 
model(^ nematicity originates from the compe- 
tition between two different stripe orders. An- 
other alternative explanation was given by or- 



bital ordering. "^^ In this scenario, the interac- 
tions between orbitals drive local ordering in 
orbital occupancy, and resulting in rotational 
symmetry breaking. However, all of these con- 
sequences are similar to our approach. It is the 
FS pocket in X and Y that are renormalized 
unequally in the ground state. 

In our approach, the nematic order is in- 
terpreted using a quasi-particle picture. Be- 
cause of the instability in the direct channel, 
the FS fluctuations (the excitations of the quasi- 
particles) can yield lower energies than its orig- 
inal isotropic FS configuration. Therefore the 
system tends to deform an isotropic FS to lower 
the ground state energy. In addition, only the 
Gaussian terms of the magnetic fluctuation were 
kept in our study. Hence, the competition with 
stripe orders does not play an essential role in 
this scenario. However, it would be interest- 
ing to study the interplay between stripes and 
nematic order (A~) by including fluctuations 
beyond the Gaussian approximation. Further- 
more, the work in this paper is based on an itin- 
erant picture. How the local ordering relates to 
the approach of this work is another interesting 
topic to be explored in the future. 



C. Concluding remarks 

In summary, we have discussed possible in- 
stabilities of the isotropic metallic phase within 
an itinerant model. By using the Hubbard- 
Stratonovich transformation, several auxiliary 
bosonic fields that describe the magnetic and 
nematic order collective modes were introduced. 
By perturbing the action of these bosonic fields 
in the isotropic phase, if the inter electron 
pocket interaction (ug) is greater than the in- 
tra electron pocket interaction (U4 ), the mag- 
netic fluctuations can drive an instability in the 
direct channels. From mean field analysis, we 
argue that nematic order is a favored ground 
state which spontaneously breaks C4 symmetry 
but preserves translational symmetry. In the 
nematic ordered state, the size of the FS pock- 
ets in X and Y becomes unequal. For A^ > 
(A~ < 0), the X (Y) pocket is larger than the 
Y (X) pocket. Because of this distortion, the 
magnetic fluctuations also exhibit anisotropy in 



X and Y directions. A similar nematic order 
can possibly occur in FS hole pocket, if angular 
depend in the coupling constants is considered. 
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Appendix A: Hubburb-Stratonovich 
Transformation 



The identity of the 'fat' unity of the Gaussian 
functional integral of 0-field is given by 



P[0] exp 



^E 



-''spm 



-M^ • M^, 



ivf— A+A+ +— A-A- 



(Al) 

Here, Un must be positive to ensure the integral 
in (IA1[) is convergent, and the 1 ('fat' unity) 
represents the normalization constant. 

The quartic fermonic interaction terms in 
(jAll) can be generated by shifting the auxiliary 
fields by 



(A2) 
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a 1 "s 
1 + 4 
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[pf 
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'g 
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wher_e s^ = Efe V^fe+,.s'^«^'V'fc,': P% = 
EfcV'fc+5,sV'L> and k = (iw„,k). w„ is the 
Fermionic Matsubara frequency. Since mJ^ > 0, 
the shift in A+ includes an extra i factor. This 
is ensures that the generated quartic fermonic 
interaction terms carry opposite signs and can- 
cel with those in ^. 
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